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Abstract 

In this letter, expectation values of exponential fields in the 2-dimensional Euclidean sine-Gordon 
field theory are calculated with variational perturbation approach up to the second order. Our 
numerical analysis indicates that for not large values of the exponential-field parameter a, our 
results agree very well with the exact formula conjectured by Lukyanov and Zamolodchikov in 
Nucl. Phys. B 493, 571 (1997). 
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I. INTRODUCTION 



This letter briefly reports our investigation on vacuum expectation value (VEV) of the ex- 
ponential field e ia ^ in the sine-Gordon (sG) field theory, G a , with variational perturbation 
theory (VPT). Here, a is a parameter, and <p(r) is the field operator at the 2-dimensional 
Euclidean space (2DES) point r = (x,t) (r is the Euclidean time). 

In 1997, starting from the exact expressions for the three special cases: the coupling (3 — ► 
(semi-classical limit), (3 = ~ and a = (3, Lukyanov and Zamolodchikov guessed an exact 
formula for the VEV of e ia ^ in the sG field theory at any (3 2 < 1 and \Re{a)\ < 1/(2/3) Q. 
Then, defining "fully connected" one-point functions, am (n is any natural number), from 
the VEV's of even-power fields <p 2n , they showed that 02 and 04 from the above-mentioned 
exact formula agree with those from perturbation theory for the sG field theory up to (3 A 
and that 02 agrees with the corresponding one-point function from perturbation theory up 
to g the coupling in the massive Thirring model, which is the fermion version of the sG field 
theory Q] . Furthermore, based on the reflection relations with Liouville reflection amplitude 
[sj], some extra arguments for supporting the conjectured exact formula were presented by 
their collaborators with them in the subsequent papers J|. Slightly later, in 2000, checks 
from perturbation theories in both an angular and a radial quantization approaches for the 
massive Thirring model indicated that the perturbation result up to g exactly coincides 
with the corresponding result obtained by expanding the exact formula according to the 
coupling g & Additional a numerical study foe the mode, in finite vo lum e also provides 
evidence for the the conjectured exact formula pj. Thus, these investigations have given 
the conjectured exact formula for the VEV of a_ complete check for the case of j3 — > \ 



(t; — ^ is equivalent to (3 2 — > |, from Eq.(8) in Ref. [l|), and some indirect evidences for its 
validity. Obviously, a direct check for the cases of (3 2 7^ ~ (Ref. [j| has provided a partial 
check for the case of j3 — > 0) is still needed. 

In fact, VPT |3] may provide such a check. It is some kind of expansion theory similar 
to the perturbation theory. However, because it properly "grafts" the variational approach 
onto the perturbation theory, the VPT produces non-perturbative results which are valid 
for any coupling strength, including both weak and strong couplings. Moreover, because 
the principle of minimal sensitivity (PMS) [8j is used to determine a man-made parameter, 
the VPT is believed to be a convergent theory, and its approximate accuracy can be sys- 
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tematically controlled and improved to tend towards the exact Q|. The primitive idea of 
VPT is not a new one and can date back to 1955 at least [9J. Now, it has developed with 
many equivalent practical schemes, and has been applied to quantum field theory (QFT), 
condensed matter physics, statical mechanics, chemical phys.es, and so on & In this 
letter, based on the variational perturbation scheme in one of our former joint pa per s [10] 
(the scheme in Ref. 10] was stemmed from the Okopinska's optimized expansion [ljj, and 



proposed by Stancu and Stevenson |l2j), we will develop a variational perturbation scheme 
for the purpose of the present letter and calculate approximately G a up to the second or- 
der. No explicit divergences exist in the resultant expression o wing to the adoption of the 
Coleman's normal-ordering renormalization prescription One will see that our 

investigation strongly support the conjectured exact formula. 

Next section, we will develop the VPT to calculate the VEV of e ia<t>{ ^ and give the VEV 
of e m ^ r ^ up to the second order in VPT. In section III, we will report our numerical results 
and make a numerical comparison between our results and the exact ones. A brief conclusion 
will be made in Sect. V. 



II. VPT FOR VEV'S AND APPROXIMATE G a UP TO THE SECOND ORDER 

We consider the 2-dimensional Euclidean sG field theory with the following Lagrangian 
density 

C-sG — 2^ V0rV V0r — 2fi COs(v^87T/50^) . (1) 

In this letter, the subscript r represents the coordinate argument, for example, = 0(r), 
and Vr 1 is the gradient in the 2DES. The Lagrangian density Eq.(l) is nothing but Eq.(5) 
in Ref. jl| if one makes the transform — > -4= (hereafter, we will use e %y fi™ a <i>( r ) instead 
of e ia< ^ and consequently the parameter a and the coupling (3 in this letter are identical 
to those in Ref. [jj, respectively). If taking \/8nj3 — > (3 and 2Q = m 2 //3 2 and adding 
the term m 2 /f3 2 in the Lagrangian density, one will get the Euclidean version of the sG 



Lagrangian density which discussed in Ref. 10]. In Eq.(l), the dimensionless (3 is the 
coupling parameter and Q is another parameter with the dimension [lenth] -2 in the natural 
unit system. It is always viable to have (3 > without loss of generality. The classical 
potential V(<p?) = — 2VL cos(v^87r/50^) is invariant under the transform — > <fi + ^|=^ with 
any integer n, and so the classical vacua are infinitely degenerate. So do the quantum vacua 



3 



of the sG field theory according to Ref . fl(| . Here, as did in Ref. [lj], we choose to consider the 
symmetry vacuum with the expectation value of the sG field operator <p r vanishing instead 
of those spontaneous symmetry broken vacua. 

The VEV of the exponential field e iV ^ a ^ is defined as follows Q 



r = . >/ftF«*(o) . = / v< t> exp{2y / 8^a0(O)} exp{- / d 2 rC sG } 

For simplicity, the exponential field in Eq.(2) is taken at r = 0. It is evident that the 
numerator and denominator in the right hand side of Eq.(2) can easily be got from the 
following sG generating functional 

Z[J] = J X?0exp{- J d 2 r [C sG - Jtf?]} , (3) 

by taking J? = i\f&Ka5(r) (5(f) = 5(x)5(r)) and J = 0, respectively. In Eq.f3), J? is an 



external source at f. For renormalizing G a , we will use its normal-ordering form Il4|. 

To perform a variational perturbation expansion on G a , now we modify Z[J] in Eq.(3) 
by following Ref. only without shifting the field </v (This is not necessary here because 
we choose to consider the symmetrical vacuum as aforementioned). That is, first introduce 
a parameter \i by adding a vanishing term J d 2 r\(j>?(ii 2 — [i 2 )4>r into the exponent of the 
functional integral in Eq.(3), then rearrange the exponent into a free-field part plus a new 
interacting part, and finally insert a formal expansion factor e before the interacting part. 
Consequently, one has 

Z[J]^Z[J;e] = exp{- J d 2 f [\l m (^ 2 ) - \l (0) (M 2 ) + \m 2 I (1) (M 2 )]} 

xexp{-e J dVW / (^-,/i)}exp{ijr 1 J} (4) 



with 



HMr, fi) = -\^ 2 <Pl - 2tt cos(v / 8^/50 r -) exp{47r/3 2 / (1) (A^ 2 )} . (5) 



2 

Here, M. is a normal-ordering mass, 



)2 , = J / (2tt)2 (p 2 +Q 2 )™ ' f° r 71 ^ 

J^ln(p 2 + Q 2 ), for n = 



hn)(Q 2 ) = 

with p a Euclidean momentum in 2DES, and J f~ 1 J = J d 2 r'd 2 r" Jff^pi J?» with 



In Eq.(6), K n (z) is the nth-order modified Bessel function of the second kind. Note that 
Eq.(4) in the extrapolating case of e = 1 is only the normal-ordering expression of Eq.(3). 

Expanding exp{— e J cPr 7~ti{jjz, A*)} i n Eq.(4) with Taylor series of the exponential, one 
can write normal-ordered G a as 



Q a _ e ina 2 I w (M 2 ) 



oo ,n(~l)" 



=0 ' 



/ nLi d 2 Wi{j§-, «) exp{I J/- 1 J}] 



87ra5(r) 



I J=0 



• (7) 



6=1 



En=o ^ / nLi ^ 2 W/(^T, A*) exp{| J/- 1 J}] 

' k 

Thus, according to the formula 0.313 on page 14 in Ref. [3], G a can be expanded as the 
following series of e 



G a = [Gf + e^W + e 2 Gi 2 ) + • • • + ^G?) + • • • ] £=1 ■ 



(8) 



This series with e = 1 is independent of the man-made parameter «, but to get its closed 
form is beyond our ability. Hence one can truncate it at any order of e to approximate it and 
then the truncated results will be dependent upon u. This arbitrary parameter u should 

n 

be determined according to the PMS |8[ as mentioned in the introduction. That is, under 
the PMS, u will be chosen from roots which make the first (or second) derivative of the 
truncated result with respect to u vanish P, 3, Q . Thus, « will depend on the truncated 
order. It is believed that it is this dependence that makes the truncated result approach the 
exact G a order by order i7|. Thus the above procedure provides an approximate method of 
calculating G a which can systematically control its approximation accuracy. It is evident 
that this procedure has no limits to the model coupling and is a non-perturbative method. 

Executing the above procedure, we have truncated the series in Eq.(8) at the second 
order of e, and the three coefficients are 



G a 0) =exp{4na 2 I {1) (M 2 )} 



expi-Jf-'J} 



(9) 



Jft=iV8TraS(r) 



exp{47ra 2 J (1) (X 2 )} 



J dV-W / (^-,//)exp{ijr 1 J} 



Jp=iV8TraS (r ) 



+G (0) 



(10) 



J~=0 
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and 



G (2) _ 1 W/ (1) (A4 2 ) 

a " 2! 



J^=i\/87ra5(r) 



2! a 



/ 

J d 2 rH I (^ J -,p)exp{^Jr 1 J} 



(11) 



J P =0 



respectively. 

Performing carefully operations in Eqs.(9),(10) and (11), we obtain the following expres- 
sion of G a approximated up to the second order of e, G T J = [G^ + eG^ + e 2 G^] e =i, 



G 



n 



V 2 W 



where, 



M 2 ' 
1 



(1 - 2a 2 K 02 ) + 4tt -^(^f^ 2 - 1 ^ 



-a 



A* ^a 2 



M 2 M 2 ' 
K 0111 + 2a\^f(K 02 ) 2 



2tt 2 y M 2J ' X M 2 

,2 



o_ ^ n 2( ^ Na 2 +/3 2 -l^ 2 ^ o2( f 1 x a 2+^_ 1 



02c 



1 " aPi^- f^Kons + Sn 2 (^) 2 (^-r 2+2 ^- 2 (K 0c ) 2 



7iM 2 ' y M 2 ' X M 2 ' y M 2 

o a 2 o i# 2 

, / " \2( P Na 2 +2 / 9 2 -2 7>- , / " \2( P Na 2 +2/3 2 -2 7^ 



Oeell 



K 02 = / dxxK$(x), 



K 0c = / <ia;a;[cosh(4a/5i^o( a; )) — 1] , 
'o 



(12) 



-^02c 

-f^Olls 
Koecll 



p2ix 1*271 poo poo 

I d9 1 / d9 2 / d Pl / d P2Pl p 2 K (R)K ( Pl )K ( P2 ) , 
o Jo Jo Jo 

p2ix poo poo 

d6 l / d9 2 / d Pl / rfp 2 p 1 p 2 (ir ( J R)) 2 [cosh(4a^ (p2)) - 1] , 
Jo Jo Jo 



p2i\ p2n poo poo 

/ d9i / d6 2 I dp x \ dp 2 p 1 p 2 K (R)K (p 1 )smh(4af3K (p 2 )) , 
Jo Jo Jo Jo 

/*27r /*27r /»oo /»oo 

/ d6 1 / d# 2 / / dp 2/Ol p 2 [exp{-4/3 2 ^ o 0R)}-l] 

Jo Jo Jo Jo 

[cosh(4a/3(Ko(pi) + ^o(p 2 )))-l] , 

/*2-7r /*27r /»oo poo 

/ d0x / dfl 2 / d Pl dp 2Pl p 2 [exp{Ap 2 K (R)}-l] 
Jo Jo Jo Jo 



Oeell — 



[exp{-4a/3(K (p 1 )-^o(p2))}-l] 
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with R = yjp\- 2pip 2 cos(#i - 6 2 ) + p|. 

In the right hand side of Eq.(12), the first two terms are the expression of G a approximated 
up to the first order of e, G T a = [G^ + eGa ] e =i- To determine the arbitrary parameter p for 



G* a , as stated in the above, we can require that = according to the PMS, and have 



d(p 2 ) 



M~ 2 = { A7r M^ K ° c a^l-2a^K 02 ) ) " (13) 

Thus, substituting last equation into the expression of G^ gives approximately the result 
of G a up to the first order in VPT. Note that M. can be taken as any positive value, and 
usually it can be referred to as unit when one renders various quantities dimensionless for 
numerical calculations. 

As for p at the second order, following the next section, one can numerically check that 
the condition ^jfrj = does not produce a real p 2 for a real value of a. However, for a real 
a, G a should be real because EVE's of odd-power fields in the sG field theory vanish (See 
Eq.(2)). Hence we have to resort to ^Sfip. = for determining p at this order (note that p 
enters the scheme in its squared power p 2 ). Fortunately, p can explicitly be obtained with 
a long expression from this condition, and then the approximate G a up to the second order 
in VPT can be concretely given from Eq.(12). 

In the same way, one can approximately give G a up to higher orders. Here we do not 
continue to consider it, and next section we will carry out a numerical calculation on the 
first and the second order results. 

In passing, G^ in Eq.(8) is only the sum of all nth-order connected Feynman diagrams 
consisting of the external vertices from exp{2-y/87ra0(O)} and n vertices from 7i/((/y, p). 
Hence, one can also obtain the above results Eq.(12) by borrowing Feynman diagrammatic 
technique with the propagator of Eq.(6). 



III. NUMERICAL CALCULATIONS AND COMPARISONS 

To perform numerical calculations, we take M. = 1 and all physical quantities in Eq.(12) 
and (13) are dimensionless. Simultaneously, one can check that this treatment (taking 
Ai = 1) amounts to taking the same normalization conditions as Eqs.(6) and (16) in Ref. 
Thus, our results can be compared with the exact formula. From Eq.(12) in Ref. Q], the 
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dimensionless Q in Eqs.(12) and (13) can be written as 



n 



?rr(i - p 2 



2T( 



2-2/3 2 



(14) 



with £ 



ii 2 



— 752 . In last equation, we took the soliton mass as unit to make Q, dimensionless. 

n 

For comparisons, we also do so for the conjectured exact formula Eq.(20) in Ref. (lf,^ 2100 *) 
when it is employed. 

Now we first report numerical results up to the first order. Mathematica5.0 programm 
gives that Kq2 = 0.5 and Kq c is finite for the case of a/3 < 0.426925 which is involved in the 
range of a/3 < ~ for the conjectured exact formula. Taking f3 = 0.2 and 0.5 as examples, 
we compared our results G ! a with the conjectured exact results of Ref. PJ in Figs.l and 2, 
respectively. In Fig.l and Fig. 2, the solid curves are G^, and the dashed curves are G exact . 
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FIG. 1: Comparison between the first-order (G^, solid curve) and conjectured exact (G^ xact , 
dashed curve) sG expectation values of exponential fields at = 0.2. 

In Fig.l, G a decreases from 1 to with the increase of a. This tendency exists in the case 
of (3 < 0.4 or so. When a < 1 and a 2 + /3 2 > 1, Eq.(13) will produce a complex fi 2 and one 
should use the second derivative for determining /i. For sufficiently small (3, one needn't do 
so because G a tends to zero with the increase of a before fi 2 approaches complex values. In 
Fig. 2, G a increases almost from 1 with the increase of a. This tendency appears for other 
larger values of (3. These figures indicate that for the case of a < 0.2 or so, the first-order 
results almost completely agree with the conjectured exact results, and for larger a, our 
results differ from the exact results with about ten percents or so (at most with 20 more 
percents when a approaches values with af3 = 0.426925 satisfied, see Table I). 
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FIG. 2: Similar to Fig.l but at f3 = 0.5. 

For calculating Gj/ , we need to calculate integrals in Eq.(12) which are involved in the 
zeroth order modified Bessel functions of the second kind. Noting Gegenbauer's addition 
formula for the zeroth-order modified Bessel functions of the second kind Q| K (K) = 
I (p 1 )K (p 2 ) + ^Y.n=i cos l n (°i ~ 02)]In(pi)K n (p 2 ) with pi < p 2 (I n (pi) is the rath-order 
modified Bessel functions of the first kind), one can finish those integrals by dividing the 
plane {pi,p 2 } into two parts: one part with p\ < p 2 and the other with pi > p 2 . In performing 
the calculations, because of the oscillatory property of cos[ra(#i — 9 2 )], it is enough to truncate 
the series in the Gegenbauer's addition formula at some n. Finishing those integrals for 
(3 = 0.5 with long-time calculations of Mathematica5.0 system, we obtained the results of 
Gq 1 and accordingly compared them with G ! a and the conjectured exact results in the table. 
In Table I, A 1 = ^"flT x 100 and A" = ^'"flT' X 100. This table indicates that G" 
has an about one percent smaller difference from the conjectured exact result than G ! a at 
(3 = 0.5. We also checked the cases of smaller and larger (3 (f3 2 < 0.426925), and found the 
same conclusion. 



IV. CONCLUSION 

In this letter, in order to check the conjectured exact formula in Ref. [if, we have calcu- 
lated the sG expectation Value of the exponential field e ta ^ <yT ^ up to the second order with 
variational perturbation approach. According to our numerical results in last section, for 
not large a, the exact formula conjectured by Lukyanov and Zamolodchikov in Ref. |l| is 
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TABLE I: Comparisons: G^ 1 , G{ and G e a xact for (3 = 0.5 





riexact 


G 1 


nil 


A 1 


A /J 


n m 


f) QQQQS1 


1 00001 758 


1 00001 1 87 


003661 


003090 


0.1 




1 001847531 


1 00137859 


36542 


0.3185 


2 


993954 


1 00850388 


1 00678750 


1 4639 


1.2912 


0.3 


990960 


1 02358276 


1 01961512 


3 2920 


2.8917 


0.4 


0.995987 


1.05423835 


1.04769147 


5.8486 


5.1913 


0.5 


1.020682 


1.11390679 


1.10415979 


9.1336 


8.1786 


0.6 


1.086765 


1.22971602 


1.21627955 


13.1538 


11.9174 


0.7 


1.242866 


1.46517988 


1.44784716 


17.8872 


16.4926 


0.8 


1.634342 


2.00162788 


1.99288037 


22.4730 


21.9378 



correct for all values of j3 in the range of (3 2 < 1, and for larger a (\Re(a)\ < 1/(2/?)), our 
numerical results can be believed to support the conjectured exact formula since our method 
is an approximate one. Thus, from the existed reports in Refs. 1,0,0] and our report here, 
we believe that the conjectured exact formula in Ref. Q] is completely convincible. 

As ending the present letter, we mention some interesting problems. The normal-ordering 
prescription amounts to a renormalization procedure on the mass parameter and makes 
finite for the range of a/3 < |. If introducing additional renormalization scheme on the 
coupling (3 and the exponential-field parameter a, it will be possible to obtain a finite G a for 
the range of a(3 > ~. Furthermore, it will be also interesting to calculate the sG expectation 
values of the exponential fields on the asymmetrical vacua. On the other hand, because it 
was an important progress in calculating VEV's of local fields, the Lukyanov-Zamolodchikov 
conjecture in Ref. Q has stimulated generalizations or conjectures on the VEV's of local 
fields in some field theories Q], and hence the method in the present letter and its gener- 
alization to finite temperature case can be used to check or confirm those generalizations or 



conjectures in Ref. 



171 ] . Besides, although our results up to the second order improved 



the results up to the first order, Table I indicates that the improvement is small. Since 
the Lukyanov-Zamolodchikov conjecture can be believed to be correct, comparing it with 
higher-order contributions in VPT for the sG expectation Values of the exponential fields 
will provide the first check in QFT on the convergency of VPT. 
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